Trion ground state, excited states and absorption spectrum using electron-exciton 

basis 
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We solve the Schrodinger equation for two electrons plus one hole by writing it in the electron- 
exciton basis. The main advantage of this basis is to eliminate the exciton contribution from 
the trion energy in a natural way. The interacting electron-exciton system is treated using the 
recently developed composite boson many-body formalism which allows an exact handling of electron 
exchange. We numerically solve the resulting electron-exciton Schrodinger equation, with the exciton 
levels restricted to the lowest Is, 2s and 3s states, and we derive the trion ground state energy as 
a function of the electron-to-hole mass ratio. While our results are in reasonable agreement with 
those obtained through the best variational methods using free carrier basis, this electron-exciton 
basis is mostly suitable to easily reach the bound and unbound trion excited states. Through their 
wave functions, we here calculate the optical absorption spectrum in the presence of hot carriers 
for 2D quantum wells. We find large peaks located at the exciton levels, which are attributed to 
electron-exciton (unbound) scattering states, and small peaks identified with trion bound states. 
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I. INTRODUCTION 

Trions, i.e., charged composite particles made of two 
electrons plus one hole (X~), or two holes plus one elec- 
tron (X + ), were first predicted in bulk semiconductors a 
few decades agc-k These 3-body objects have well-known 
counterparts in atomic physics: charged hydrogen ion 2 - — 
H~ or molecule H ^ ■ Although their binding energies are 
several orders of magnitude larger than for semiconduc- 
tor trions, their physics are essentially the same provided 
that we use semiconductor effective masses, i.e., we ne- 
glect all complexities coming from the valence band by 
assuming one hole species only. 

In the late 70's, trions with very small binding en- 
ergy has been observed in bulk semiconductor sam- 
ples such as Ge and Si (Refs. |6||7|). However, it is in 
the 90 's only that high-quality semiconductor quantum 
wells allowed precise studies of trion states using op- 
tical measurements^—. Indeed, the reduction of di- 
mensionality is known to enhance all composite particle 
binding energies. Recent experiments in highly doped 
materials^ 2 - - — have further investigated the interesting 
interaction which exists between trions and free carriers. 

Although single trion eigenstates have up to now been 
approached through a free carrier basis made of two elec- 
trons plus one hol o 15 i 16 — or two holes plus one electron — 
this basis is not the most appropriate one from a physical 
point of view. Indeed, the three free carriers of a trion 
interact through the strong Coulomb interaction which 
exists between elementary charges ±e. This interaction 
has to be treated exactly in order to possibly find the 
poles associated with the trion bound states. It is phys- 
ically clear that the trion ground state energy must lie 
slightly below the exciton energy, because the "second" 
electron is attracted not by an elementary charge but by 



the excitonic dipole. Yet, this weak attraction still has to 
be treated exactly in order to get the trion bound states. 
Thus, in the approach using free carriers, we face the 
challenging task of solving a 3-body Schrodinger equa- 
tion with a very high precision, since ultimately the trion 
binding energy we arc interested in, only is a very small 
fraction of the total trion energy we get — the trion bind- 
ing energy being the difference between the trion and the 
exciton ground state energies. 

The procedure we here propose avoids this difficulty. It 
uses a basis made of electron-exciton pairs (e, X), the ex- 
citon states being analytically known, both in 2D and 3D. 
Instead of a free energy made of free carrier kinetic ener- 



gies 



e k^ W1 th a minimum value equal to zero, 



the (e, X) basis we here use already contains the exciton 
binding energy at zeroth order in the interaction, since 
the "free energy" then is Exq+s^) ■ The construction of 
a trion Schrodinger equation using electron-exciton pairs 
as a basis can be done by using the recently developed 
composite boson many-body theory^ that can exactly 
handle carrier exchanges between exciton and free car- 
rier as induced by the Pauli exclusion principle. 

The goal of this paper is not to compete with the vari- 
ational approaches commonly used to derive the trion 
ground state energy. Thanks to rapid advance in com- 
puter technology, elaborate variational methods can now 
reach the ground state energy with an amazing precision 
of ten decimal digits or morei^. Rather, our purpose in 
proposing a new way to approach the trion problem is 
threefold: 

a) Variational methods are quite appropriate to derive 
the trion ground state energy with a numerical accuracy 
which surpasses experimental data. However, the coor- 
dinates commonly usecU^ in these variational procedure 
completely mask the underlying physics of an electron 
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interacting with an exciton through Coulomb processes 
and the effect of carrier exchanges on the structure. In 
contrast, the method presented here provides a physi- 
cally clear picture of the many-body effects involved in a 
trion. Moreover, we can reach a reasonably good value 
of the trion ground state energy by numerically solving 
the trion Schrddinger equation with just a few low-lying 
exciton levels kept in the electron-exciton basis, which 
makes this numerical resolution quite easy to perform. 

b) The trion operator, defined in Eq. (fT0|) . constructed 
on the electron-exciton basis we here use, is quite suitable 
for further study of interacting trion a 20 ' 21 because it inte- 
grates the trion into the general framework of interacting 
excitons, for which a powerful many-body formalism has 
been developed^. Solving the trion Schrddinger equa- 
tion in the electron-exciton basis, as addressed in the 
present work, provides useful grounds for future works 
on many-body effects involving trions. 

(c) Finally, the electron-exciton basis allows us to ob- 
tain the trion excited states as easily as the ground state, 
these excited states being difficult to reach within varia- 
tional procedures. 

Using these ground and excited state wave functions 
makes it possible to compute the whole absorption spec- 
trum of bound and unbound trions in doped semicon- 
ductors. We here concentrate on optical absorption in 
2D quantum wells in the presence of hot carriers. For 
X~ trion, i.e., two electrons and one hole with electron 
masses smaller than the hole mass, this spectrum shows 
one small peak below the ground state exciton energy, 
which is associated with the trion ground state, and right 
at the Is, 2s and 3s exciton levels, larger peaks associ- 
ated with unbound trions, i.e., electron-exciton scatter- 
ing states. Through the carrier energy distribution, the 
absorption spectrum actually acts as a tool to probe the 
trion relative motion wave function in momentum space. 
The peak height is found to decrease with photon energy, 
due to the fact that the trion wave function amplitude 
decreases when the electron relative motion momentum 
increases (see Figs. ITT1 and H"2"|) . 

The paper is organized as follows: 

In Sec|TTl we write the trion Schrddinger equation using 
the two-electron-one-hole basis for completeness, since 
this basis is the standard one to approach trion. We 
pay particular attention to the parity condition for trions 
made of two electrons with same spin or opposite spins. 

In Sec lIIIl we write the trion states in the electron- 
exciton basis and we show how the parity condition ap- 
pears within this (e, X) representation. We then derive 
the Schrddinger equation fulfilled by the prefactor of the 
trion expansion in this (e, X) basis. In the last part, we 
focus on trion made of opposite spin electrons in a sin- 
glet state to possibly reach the trion ground state and 
also the trion relative motion wave functions appearing 
in photon absorption spectrum. 

In Sec lIVl we first present numerical results on the 
binding energies of the trion ground and excited states for 
various electron-to-hole mass ratios. We then show the 



trion relative motion wave functions for bound and un- 
bound states. Finally, we use these trion wave functions 
to calculate the optical absorption spectra for various hot 
carrier distributions. 
We then conclude. 



II. X~ TRION IN TERMS OF TWO 
ELECTRONS PLUS ONE HOLE (e,e',h) 

We consider two electrons with spins s and s' and one 
hole with total angular momentum m. The electron spins 
can be ±1/2. In bulk samples, the hole angular momen- 
tum can be m = (±3/2, ±1/2), while in narrow quantum 
wells, it reduces to m = ±3/2 due to the heavy- light hole 
energy splitting induced by the well confinement. How- 
ever, since the hole angular momentum does not play a 
role for a single X~ trion if we neglect the exciton split- 
ting due to electron-hole exchange, we will in the follow- 
ing ignore the hole angular momentum m for simplicity. 

The basis commonly used to represent X~ trion is 
made of three free carrier states, namely ^a^ 
By noting that the system Hamiltonian for two electrons 
(e, e') plus one hole (h) in first quantization, 
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(1) 



is invariant under the (e O e') exchange, where e sc is 
the semiconductor dielectric constant. We conclude that 
the orbital part of the trion wave function must be even 
or odd with respect to this exchange — the ground state 
wave function being even in order for this state to be 
nondegencrate. As a result, the trion eigenstates can be 
divided into two groups. One group consists of states 
that have an even orbital part and odd spin part (due 
to the Pauli exclusion principle). Those are spin singlet 
\S = 0, S z = 0). The other group consists of states that 
have an odd orbital part and even spin part. Those are 
spin triplets \S = 1,S Z = (0,±1)). 

Consequently, the eigenstates of the Schrddinger equa- 
tion for one trion in the (e, e', h) basis, 



(2) 



can be written as 



(T) 

S=l,S z =2s 



E 

k c k'k h 



^kek'k. 



(3) 



when the two electrons have the same spin s, and as 



IT) 



u k c ;l/2"k; ; -l/2 



' — E ^k e k^k h 
kek^k,, 

-(-l) S <;-l/24 i; l/ 2 | bljv) (4) 
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when the two electron spins are opposite. 



Since a k . s a k , , = —a\j .y^k . 8 , it is possible to replace 

the prefactors in Eq. © by (i&Sg^ - Wfkfk^) /2- 
For the same reason, we can rewrite the second operator 
in Eq. (j4]) as the first one, the prefactors now reading as 
(— l) S Vvkfk ■ As a resu lt, instead of writing the four 
trion states as in Eqs. (|3I4[) , we may as well write them 
in a more compact form as 



(T) 

S,S z =s+s< 



= E €AA.A>' b U^ (5) 



k e k'k h 



where the prefactors in this expansion now read 



6 (T,5) 
u k e k'k,, 



within an irrelevant normalization factor. This new pref- 
actor then fulfills the parity condition 



&l = (-i)**SfL 



(7) 



a requirement not necessary to enforce if we use the ex- 
pansions pi4[) with prefactors V'k k/ k ' smce this parity 
condition is then fulfilled by the operators themselves. 

Due to the Schrodinger equation ([2]), the prefactors of 
the trion state in Eq. (|5|) must fulfill 



= 



E {K^-£ (r,s) 

^E^Kr+ik- q , kh -^ 



J k c k'k h 



(T,S) 

k^+q.k^k^-q 



b {T,S) 

k e ,k;,+q,k h -q 



(8) 



where the free part of the energy reads as -Ek e k' k h 



e (e) + e (e) 

fc k c + fc k' 



If we now project this equation onto 



{v\b Ph a v ' e . s >a v<i . a and use 

(«|ap'e;s' a Pe;s«k e ;s a k^s'l w ) = V^'^Pek 

we get one term only for s' = 
(0, 1), S z — 0) trions, namely 



<5ss"5p^k e ^p c k^ 

s, i.e., for (5 = 







,(T,S) 

Pe+q.p'e-q.Ph 



6 (T,S) 

Pe+q,p' e ,Ph-q 



^PciPe+qiPh-q 



(9) 



By contrast, for s' = s, i.e., for (S = 1,S Z = 2s) triplet 
trions, this projection yields two terms. These two terms, 

{T S— 1) [T S— 1) 

however, are equal, since k ,.k'k h = — ^k'k^kh arl( l 
Vq = VL q , so that we end up with the same Schrodinger 
equation © for all four trion states. 



It is actually possible to recover this Schrodinger equa- 
tion starting from Eqs. (|3I4|) in which the trion states 

are written in terms of ip^ k , kh , instead of using Eq. ([5]) 

(T S) 

with prefactors k k , k . The derivation can be found in 
|Appeiidix A 



The energies of trion states using the (e,e',h) basis 
then follow from solving a 3-body Schrodinger equation 

(HJ) with the parity condition 4>k'^u h = (-l^vkjc,,- 
The trion ground state belongs to the set of even states, 
i.e., states with S = 0, while the solutions with odd wave 
functions lead to triplet states, these three states being 
degenerate. 



III. TRION IN TERMS OF 
ELECTRO N-EXCITON PAIRS (e,X) 

We now turn to the trion representation in terms of 
electron-exciton pairs. Some useful results on excitons 
interacting with free electrons are briefly rederived in 
| Appendix B| References to equations in this Appendix 
start with the letter B. 

Since excitons form a complete basis for (e, h) pairs, 
we can rewrite the trion states of Eq. ([5]) as 



(T) 

S,S z =s+s' 



E 

k, i 



Bl 



t k e ;s- 1 -'i;s' 



(10) 



where Eq. (|B.2j) allows us to relate the prefactors of this 
(e, X) expansion to the ones of the (e, e', h) expansion in 
Eq.© through 



k' c k h 



(11) 



Before going further, let us note that the parity con- 
dition ([7]) on </> k k f k , resulting from the Pauli exclusion 
principle between the two electrons of the trion, imposes 



E A (5 



k, i 



EE E (i'Kci h )(ci h K\i)(i\PePh)^JL 

k c i qh PcPh 

E Ck?L<*'l k eP/>} = ("I)* E *S?p fc <i'|keP* 
k e p h 



kePh 



(12) 



where A( ? ' I is the Pauli scattering for electron ex- 

v i i ' 

change between a k e electron and a i exciton, defined in 
Eq. (|B.4[) and shown in the diagram of Fig. [1] So, the 
parity condition ([7]) reads in the (e, X) subspace as 



$Zf = (-d s E^ 



(13) 



k, i 



4 




/ k k \ 

FIG. 1: Pauli scattering A( ., e . J for electron exchange be- 
tween a k e electron and a i exciton, as given in Eq. (IB.4II . 



gives, using Eq. (|B.9 



o=Ki-^ (r - s) )C;f 



. (i6) 



If instead, we project Eq. ([15)) onto («|Bj ; i/2a Pt ,;-i/2) we 
get a somewhat more complicated equation, namely, 



o = X>("; k ; 



k e ,i 



EfTrfjigf (17) 



which also reads as 



k\ 




i = {y i ,Q i ) 



-P + a h Qt 



/ k k \ 

FIG. 2: Pauli scattering A( e ) for electron exchange be- 
tween a k e electron and a i = (yi, Qi) exciton, as given in 
Eq. MM - 



which clearly is more complicated than Eq. ([7]). We will 
show below how to overcome this complication. 

For H being the Hamiltonian of a system with ar- 
bitrary numbers of electrons and holes, we get using 
Eqs. (|R10| and (|B~Til 



HBl s ,a{jv) 



B ls>H)<Jv) 



= E^ate-XM + E 1? k ; >k-A;^m 

k' e i' 

where £k e ,i = + -^1 is the free electron-exciton pair 
energy. The electron-exciton scattering £ ( k f ke ) associ- 
ated with their direct Coulomb interactions is shown in 
the diagram of Fig. EH 

From this equation, it is easy to show that the 



I* 



(T) 

S,Sz=s+s 



i ) state in Eq. (fT0|) is eigenstate of the Hamil- 



tonian H with energy E^'^ provided that 



E 



E^ 



k'i' 



°L.4»- (is) 



(i) Let us first consider s = — s' = 1/2. The projec- 
tion of the above equation onto (i>|-Bj ; -i/2%> e ;i/2 readily 



= 



E 



Pe k e ' 



Ew e ,i — E 



(T,S) 



4 (T,S) 
k e ,i > 

(18) 

where £ in ( F ^ ek i c ) is the exchange-Coulomb scattering 
defined in Eq. (|F3. 14|) . where "in" refers to the fact 
that Coulomb interactions take place between the "in" 
electron-exciton pair (k e ,z), as shown in the diagram of 
Fig- El It is actually possible to show that Eq. (TTS)) follows 
from Eq. by using the parity condition (|13[) . Indeed, 
Eq. (fT^|) allows us to transform the RHS of Eq. (Til))) into 



(-1) 



S E 

k c i 



(is PcJ -^))A(P; k ;)+e-(P; k ; ) 



where £° 



(19) 
Using 

Eq. (|B~T8|) for the link between £ in and £ out , we then 

find that the <fii<!~'P prefactors in Eqs. (JTHJ) and (fT9"]) are 
indeed equal. 

(ii) We now consider Eq. ([T5|) for s = s' and project 



is defined in Eq. (|B.15j) . 



it onto (i>|ifj;«ap a ;g- According to Eq. (|B.9|) . we get two 
terms. These two terms just correspond to the RHS of 
Eqs. (jTSJ) and (p~7|) . so that we again end up with <K, e 'j 



fulfilling Eq. ([16)) for S z = ±1. This is fully reasonable 
because Eq. (fTB)) holds for triplet and singlet trions made 
of opposite spin electrons, while the three triplet states 
(S = 1,S Z = (0, ±1)) are degenerate. So, they must obey 
the same Schrodingcr equation. 

As a result, the Schrodingcr equation for one trion in 
the electron-exciton basis reads for all trion states as 







k e ,i 



E (T,S) 



k'i' 



k e k' 



(20) 



where the 



wave function in this basis is subject to 



the constraint imposed by the parity condition () 13[) - The 
singlet and triplet trion states are given by Eq. (TTU)) with 
the prefactor (j)^ '. fulfilling the Schrodingcr equation 

The parity condition (|13[) . however, is difficult to im- 
plement numerically. It is actually possible to include 
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>K-K 

Pe . 



X- 




kl 



a: 



/ k' k \ 

FIG. 3: Direct Coulomb scattering H " ) between a k e 

° ^ V i . i L. 

electron and a i exciton, as given in Eq. (|B.12|) . 



/ k' k \ 

FIG. 5: "In" exchange- Coulomb scattering f ln ( , e ) b e- 
tween a k e electron and a i exciton, as given in Eq. (|B.16I) . 



k'., 



P+aMi \ P + a eQi 

i' = (v,,a)_< • _<_ i = (v,,a) 



k' 




■ P + <* h Qe 



f k' k \ 

FIG. 4: Part of the direct Coulomb scattering £( 6 J 
between a k c electron and a i = (^i,Qi) exciton, as 
given in Eq. (|B.13[) . coming from electron-electron repulsion. 

/ k' k \ 

£( . ) also contains a negative contribution coming from 
electron-hole attraction, as seen from the second diagram of 
Fig.H 



this condition into the Schrodinger equation itself by set- 
ting 



dT,S) _ (T,S) 



^T,f + (-i) s E A ( k ; ■' 1 )^ r,s) - 



(21) 



k'i' 



so, the parity condition (|13[) on <fiZ'P is automatically 
satisfied for arbitrary ip^ '■ , which can then be re- 
garded as a free function. By inserting Eq. (|21[) into the 
Schrodinger equation (|2T)|) . we find the equation fulfilled 
by the fZ'i function as 



FIG. 6: Part of the "in" exchange-Coulomb scattering 

k' k 

£ m ( i'" ■ ) between a k e electron and a i = (vi,Qi) exci- 
ton, as given in Eq. (|B.17|) . coming from electron-electron 
repulsion. 



where the relevant set of exchange processes now take 
a symmetrical form between the "in" state (k' e ,i') and 
"out" state (k e ,i), namely 
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cch \ 



k e k' 



k c k'_ 



out / k e k' 



E K ,i>)X 



k e k' 



(24) 



k'i' 



k'i' 



(22) 



Using Eq. (|B.18|) . it is possible to rewrite this equation 
into a symmetrical form with respect to the (k e ,i) and 
(k' e , i') states, as 



e^ZT +£[C e ^) + (- 1 ) s ^ 



We note that this effective exchange scattering not only 
contains the "in" and "out" Coulomb-exchange scatter- 
ings, but also the Pauli scattering multiplied by the sum 
of "in" and "out" electron-exciton pair energies in order 
to end with an energy-like quantity. 

From now on, we only consider trions with zero center- 
of-mass momentum. These arc constructed on zero- 
momentum electron-exciton pairs, k and i = (is, — k). 
Moreover, we concentrate on singlet trion states, S = 0, 
which are the relevant states for the determination of 
both the trion ground state and the photoabsorption 
spectrum in the presence of a doped electron gas. So, 
we will not quote the S index anymore to simplify the 
notation. 

Equation (|23[) can be formally written into a matrix 
form as 



k e k^ 

i i' 



k'i' 



fZf +(-i) s E*( 

k'i' 



(T,S) 
^k'.i' 



(23) 



(25) 



where E^ is the trion energy we want to determine. 
For zero center-of-mass momentum pairs, as considered 



6 



from now on, the components of the ifP^ vector are 
^(k, v) = V^Lfiz-k)- The identity matrix is dc- 



noted as 1 while E is a diagonal matrix with components 
(kv\E\k.'v') = E(k, v)5\^i5 vu i , the free electron-exciton 
pair energy being E(k,v) = £ k ,(, y _ k ) = e v + \i 2 /2fi eX 
where = to" 1 + (m e + TOfe) -1 is the inverse of the 
electron-exciton pair reduced mass. £ e g is an energy-like 
matrix which includes Coulomb and exchange processes 
within the electron-exciton pair. Its components are 

(H&fflkV) = [e( (!/ . k k) (,' k - k '))+^ h ( ( ,. k k) (l / k - k <) 

(26) 

Finally, A is a dimensionless matrix which originates from 
the fact that the hole in the electron-hole pair forming the 
exciton can be associated with any of the two electrons. 
Its components are (ki/|A|kV) = A(^ k k ^ ( I / k _ k '))- 

Figures 1-4 show the various scatterings appearing be- 
tween an electron and an exciton (for their derivations, 
see |Appcndix B I . For zero-momentum electron-exciton 
pairs, they take a compact form: The electron exchange 
reduces to 




0.05 



FIG. 7: (Color online) Ground state trion binding energy in 



(2D) 



unit of the 2D exciton binding energy R x 
a function of the electron-to-hole mass ratio 
found for m e slightly lighter than mn- 



4i? 



(3D) 



a.s 



A minimum is 



M(, k 



(u-k) (i/.-k') 



= (z/|k' + a e k)(k + a e k'| 



(27) 



while the direct Coulomb scattering between an electron 
and an exciton simply reads as 



k' 



£( K-k) (i/',-k') . 



(28) 



Vk- k < J2 [HP - «/»k)(P - «fckV> 
p 



(oth 



When Coulomb scattering is mixed with electron ex- 
change, the resulting exchange- Coulomb scattering is 
given, for zero-momentum pairs, by 




k k' 

(i>,-k) (v',-k') , 



e 



k' k 

. {v',-k') (i/,-k) . 



^^k'+a e k-p(^|p) 
P 

x[(p + a/l (k-k>') - (k- 



,k'| 



(29) 



where Vq — 2ire 2 / t sc L 2 q in 2D and = 4ne 2 /e sc L 3 q 2 
in 3D. 

It is actually possible to cast Eq. (1231) into a 
Schrodinger-likc equation for ip^\ by rewriting it as 

£(ryr) = (i + X)-\e + i s )v {r) = iW r) - (30) 

The operator H c r then appears as the effective Hamil- 
tonian for trion written in the electron-exciton basis. 
We see that electron exchange within this pair appears 
not only through the energy-like exchange scatterings 
(£ m ,£ out ) and (EX) included in £ ec h, but also in a more 
subtle way through the (1 + A) -1 operator appearing in 
front of the more naive part of the Hamiltonian, namely 
(E + is). 



FIG. 8: (Color online) Same as Fig.[7]for 3D trion, the binding 
energy unit being R^ D ^ . The behavior is qualitatively similar, 
but the binding energies are significantly smaller, which makes 
the experimental observation of 3D trions more difficult. 



IV. RESULTS AND DISCUSSIONS 



The last step is to numerically solve the Schrddingcr 
equation (|2"5|l or cquivalcntly Eq. (|30|) . To do so, we 
restrict the exciton basis (p|f) to the three low- lying s- 
like states. The inclusion of p- or d-like exciton wave 
functions requires a very careful treatment of the angular 
part of the electron-exciton pair wave functions in the 
calculation of the A and £ scatterings, which is beyond 
the scope of the present work. This extension will be 
considered elsewhere. The 2D normalized wave functions 
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(p\v) used in this electron-exciton basis are22 



(f 



'2tt 



(1 + a 2 f p 2 /4) 3 / 2 ' 



<p|2s) = 



/ax\ 3VfoF(9n^p 2 /4 - 1) 



U/ (l + Oa^p 2 /^ 5 / 



(31) 



ra x \ 5 VlMl - 25(4 p 2 + 625a^p 4 /16) 



(P|3S> \L) (1 + 25a^p 2 /4)V2 

while for 3D they read as 



<P|1*> 
(P|2s) 
(p|3a) 



1) 



/ax_\ 3 / 2 80r 
\LJ (l + a^p 2 ) 2 ' 
/a x \3/2 32 % /2T(4a 2 t p 2 
\T) (1 + 4a 2 . p 2 

/a x \3/2 72y / 37r(l - 30a 
VIT/ (TTiTa^p 



2 2 



81a 4 , p 4 ) 



2 "\ 1 



L is the sample size, ax is the 3D exciton Bohr radius 
defined as ax = S 2 e sc //ijife 2 , with \i x l = m^ 1 + yn^ 1 
being the inverse of the exciton reduced mass and e sc 
being the semiconductor dielectric constant, one order 
of magnitude larger in semiconductor samples than in 
vacuum. 

The Schrodinger equation (|30|) . seen as a matrix 
spanned by the k momentum, can be solved as a gener- 
alized eigenvalue problem. To this end, we sampled the 
k = |k| value with 100 mesh points, according to fc, = uf 
where the Uj's are equally distributed, thereby allowing 
for more sampling in the small k region. However, for 
TOg/m/j = 1/50 and 1/100 in 3D, we have taken 220 and 
350 mesh points in order for the A matrices to unam- 
biguously stay positive definite. We also took an upper 
cutoff fc max of 10 in 3D but 20 in 2D (in unit of ax), since 
the 2D wave functions have a larger radial extension in 
k-sspace. 



A. Trion binding energies for ground and excited 
states 

Figures [7] and [8] show the binding energies of the 2D 
and 3D ground state trions as a function of the electron- 
to-hole mass ratio m e /irih. These two sets of results are 
expressed in terms of their 2D and 3D effective Ryd- 



bergs, namely R^ D ^ and R^ D \ with R 



and R 



(3D) 



(2D) 
V X 



AR 



(3D) 
X 



(p,x/^ c )13.6 eV, the ratio of exciton 



reduced mass to free electron mass mo being p, x = 
m e mh/mo(m e +mh)- Although the qualitative behaviors 
of the 2D and 3D curves are very similar, it is worth not- 
ing that the absolute values of the 2D binding energies 
are significantly larger than in 3D. This is mainly due 
to the fact that the Coulomb interaction decreases 
more slowly, as 1/q in 2D instead of 1/g 2 in 3D. For 
m e /mh ^ 1, the trion binding energy stays essentially 
constant as a function of mass ratio, in agreement with 
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(32) FIG. 9: (Color online) Bound state binding energies of the 2D 



trion ground state(GS) and the first 3 excited states in R^ D ^ 
unit as a function of the electron-to-hole mass ratio. 
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FIG. 10: (Color online) Same as Fig. [9] for 3D, the energy 
unit being R^ D \ 



experimental data: in most bulk semiconductor sam- 
ples like GaAs and InAs, with an effective electron mass 
smaller than the effective hole mass, the X~ trion binding 
energy docs not significantly depend on the mass ratio. 
This is also true for 2D. 

In contrast, when m e /rrih increases above 1, which 
mimics the change from X~ to X + trions, with two car- 
riers then having a heavy mass, we see that the binding 
energy first increases essentially logarithmically and then 
saturates for large mass ratios: when m e /mh 3> 1, the 
trion behaves as an hydrogen molecule- ion rlj, with a 

strong binding energy known to be 0.412 R^ D ^ in 2D 
(Ref.il) and 0.194 R ( ^ D) in 3D (Ref.i). 

The binding energies of the 2D and 3D trion ground 
states we obtain as a function of mass ratio, are 
qualitatively similar to those obtained from variational 
methods^ 4 - - — . However, the values we find only account 
for 60% — 70% of the most accurate variational results 
when m e /m/j > 1, and 50%— 60% only when m e /mh 1- 
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This suggests that the polarization of the exciton in- 
duced by the second electron, which is not accounted for 
here because our numerical calculation excludes p-like, 
d-likc, and higher angular momentum exciton wave func- 
tions, plays a significant role in the binding, and this 
is more so when two carriers have a light mass. For 
GaAs, with effective electron mass m e = 0.067too, ef- 
fective hole mass rrih = 0.34mo, and dielectric constant 
e sc = 12.5, Fig. Ogives a 2D trion binding energy equal 
to 0.06 R^ D ^ = 1.2 meV, which is about a factor of 2 
smaller than the best variational resul t 24 ' 26 i 27 . 

In spite of this discrepancy, one important advantage 
of the present method compared to usual variational ap- 
proaches is that it allows to reach trion excited states as 
easily as the ground state, while the excited states are not 
easily obtained by usual variational methods. Figures [§] 
and [10] show the various excited bound states of 2D and 
3D trions as a function of the mass ratio m e /mh- For 
m e /mh < 1, we find one bound state only, which is the 
ground state, in agreement with known results for X~ 
trions, with a hole mass larger than the electron mass, 
both in 3D (Ref. HH) and in 2D (Ref. HJ). By contrast, 
the number of trion bound states increases with m e /mh 
when this ratio gets larger than 1. A second bound state 
emerges for me/mh ~ 16 in 3D and m e /mh ~ 10 in 2D. 
And so on... These additional bound states mainly come 
from the vibrational motion of the two heavy carriers. 
The trion then is similar to a hydrogen molecule-ion 
or its deuterium isotopes HD + and D^~. Spectroscopic 
studies of their vibrational and rotational levels have 
been a central issue in molecular physics because they 
represent the simplest molecules ever theoretically22r— 
and experimentally^— studied. For very large mass ra- 
tio m e /m/j = 1000 close to the molecular limit, we find 15 
and 12 bound states in 2D and 3D systems, respectively. 




P 



(b) 

FIG. 11: (Color online) L 2 |{r, p|r?>| 2 for the trion bound 
ground state r\ — rjo (figure a) and |(r, p|??)| 2 for an unbound 
state in a semilog plot (figure b). The mass ratio is taken equal 
to mh/rrie = 5. The units for relative motion momentum p 
and exciton radial coordinate r are ax and a^ 1 , respectively. 



B. Trion wave function 

By solving the Schrodinger equation (|3TJ)) . we can also 
obtain the ip^ _jn function from which the trion 

wave function 0ki=(y O -k) w ith proper symmetry prop- 
erty in the electron-exciton basis, can be reached through 
Eq. (|2"Tj) . Indeed, since Pauli scatterings A( ke k f ) con- 

(T 5—0) 

serve momentum and spin, the ip^, \, functions which 
are required in the sum of Eq. (|2ip . also correspond to 

(T,S=0) r 

^,i'=(^',-k' e ) unctions. 

To go further, let us first relate the 4\ 'i=( v _k ) mnc_ 
tion, which we can derive through solving the Schrodinger 
equation fl30[) . to the trion relative motion wave function 
used in previous works on trion. In Ref. l37l we have 
shown (see Eq. (3.16) in this reference) that the wave 
function of a trion with center-of-mass momentum K, 
relative motion index rj and spin S = (0, 1) splits as 

(r e> i V> r h |K I T ? ,S) = (R T |K)(r,u|r7,S), (33) 



where = (m e r e + m e r e > + m^v^)/ (2m e + rrih) is 
the center-of-mass coordinate, r = r e — is the dis- 
tance between one electron and the hole, while u = 
r e i — (m e r e + m^r^) / '(m e + rrih) is the distance between 
the other electron and the center-of-mass coordinate of 
the electron-hole pair. From (r, u|?7, S), we can con- 
struct (v, p\rj, S) through a double Fourier transform (see 
Eq. (3.18) of Ref. H3) 

(p,p\v,S) = ydrdu(Hr)(p|u)(r,u|77,S). (34) 

The parity condition (r e , r e i, r/JK, r), S) = 
(— l) s (iv, r e , r/j|K, rj, S) on the (r e ,rv) coordinates 
then leads to 

{u,p\v,S) = (-l) s ^(^|p'+ ae p)(p+« e p>')(^,P'h^) 

u',p' 

(35) 

(T S) 

which exactly is the parity condition for tfi^ '. given in 
Eq. (fl~3"|) when i = (v. — k e ), i.e., when the center-of-mass 
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FIG. 12: (Color online) | (r = 0, p\rj) | 2 for three unbound trion 
states in a semilog plot, the mass ratio being rrih/m^ = 5 and 
the momentum p unit being ax ■ 



momentum of the electron-exciton pair is zero, because 
the Pauli scattering then reduces to Eq. (f27|). 

If we now consider the creation operator for the 
(K, 77, -S) trion made of opposite spin electrons, we find 
that it reads in terms of electron-exciton pairs as (see 
Eq. (3.44) of Ref. 13) 



T, 



t 



KrjS,S*=0 



t 



p# e K,l/2 »,-p+feK-l/2' 

(36) 

So, in view of Eq. (JTUJ), we are led to conclude that 
(v,p\rj,S) must be identified with 



= 6 {T ' S) 
r p,i=0,-p) 



(37) 



with T = (K = 0,77). 

Actually, the relative motion wave functions that have 
physical relevance is not (is, p\rj, S) but (r, u|t7, S) or 
(r, p|?7, S). For bound state trion, (r, u\r), S) has an ex- 
tension of the order of the exciton size ax along r and 
an extension of the order of the trion size ax along u, 
while this latter extension goes up to the sample size L 
for unbound trions. Through a Fourier transform, we can 
reach 



r,p\rj,S) = / du(p|u)(r,u|77,5), 



(38) 



which has an extension along p of the order of 1/ay for 
bound state (see Fig. 11(a) ) and of the order of 1/L for 
unbound states (see Figs. 11(b) and[T2|). 

Due to dimensional arguments, the normalized 
|(r, u 1 770 ) 1 2 which extends over ax for r and over o/p 
for u in the case of bound state, must be such that 
l~a£a£|(r = 0,u = 0|?7o)| 2 . Since (p|u) = ^/L D ' 2 , 
Eq. ([58} then gives, using again dimensional arguments, 

(r = 0,p = 0|r ?o > ^a^yr 
trion ground state with 



x a r ! so ' we en d fo r the 



|(r = 0,p = 0|77o)p 



ax 
axL 



(39) 



In the case of unbound trions, the u extension is of the 
order of L instead of a^; so, the same dimensional argu- 
ments give 

|<r = 0,p| ?? )| 2 ^ ■ (40) 

The |(r = 0,p|r?)| 2 function is peaked on various mo- 
menta p v depending on the unbound trion energies, as 
understood by noting that, in the absence of electron- 
exciton interactions, the trion wave function (r, u\r]) must 
reduce to a product of an exciton wave function and 
an electron plane wave, (r,u\rj,S) ~ </>i,(r)(e* p,7 ' u / L D / 2 ). 
From 



r,u\v) = y]{^W){u\p)(v,p\v), 



(41) 



with (u|p) = e tp u /L D / 2 , it is then easy to see that 
this requires (v,p\rj) to be a delta function d p p , this 
delta function being possibly broadened into a peak by 
electron-exciton interactions. 

We can compute (r,p\rj,S) from the (y, p\r), S) func- 
tion that we obtain through the numerical resolution of 
Eq. (|30|) , via a "Fourier transform in the exciton sense" , 
namely 



v,p\V, s } = / dr(v\r)(r,p\ri,S), 



or equivalently 



(r, pl^, S> = 5^<r|i/><i/, p|»,, S>. 



(42) 



(43) 



Note that for r = 0, the v exciton levels that survive 
in the v sum are s-likc states only. In practice, in or- 
der to compute (r,p|?y, 5* = 0) from (1/, p|r/, S = 0) = 

dr' S 7 G " > . with T = (K = 0, ??), as obtained from the 
resolution of the Schrodinger equation (J30J) along with 
the parity condition given in Eq. (|21[) . we have only 
kept the three lowest s states in the v sum of Eq. (|43|) . 
which is consistent with the fact that we have reduced 
the electron-exciton basis in this Schrodinger equation 
to these three lowest levels. 

Figure (fTTI) shows |(r, p|?7, S — 0)| 2 for the trion bound 
ground state, 77 = 770, and for an unbound state. Note 
that, due to Eq. ([BP]) , we are forced to plot the ground 
state through £ 2 |(r, p|t7 , S = 0)| 2 in order to have a 
quantity independent of sample size L. The function 
|(r,p|?7, S = 0)| 2 , from now on written as | (r, p>| 77} | 2 for 
simplicity, is concentrated around the origin, r = 0, p — 
0, for bound state trion, whereas it is peaked around 
a p^ value for unbound states. (In practice, the 1/L 
width in this figure comes from the momentum discretiza- 
tion used in the numerical resolution) . Since the relevant 
quantity for the trion oscillator strength is going to be 
|(r = 0,p|t7)| 2 , we have also given a precise plot of these 
normalized functions for three unbound states (Fig. [T2j) . 
We find that when the unbound trion energy increases, 
this function broadens due to electron-exciton interac- 
tions. 
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C. Trion absorption spectrum 

Knowing the trion wave function for bound and un- 
bound states, it becomes possible to compute the trion 
absorption spectrum as a function of the photon en- 
ergy. Trions are usually observed in the optical spectra 
of lightly doped semiconductor quantum wells. When an 
exciton is created in a semiconductor through photoab- 
sorption in the presence of an excess of free electrons or 
free holes, the photocreated exciton interacts with a free 
carrier to form a trion which can either be in a bound or 
in an unbound state, depending on the photon energy. 

Let us first consider an initial state with one free elec- 
tron with momentum and one photon with energy uj p h 
and momentum Q p h- After photon absorption, the final 
state contains two electrons and one hole, their center- 
of-mass momentum being + Qpft. Since we are mainly 
interested in the low-lying trion states, we here focus on 
the spin singlet (S = 0) states. The Fermi golden rule 
gives the photon absorption as (—2) times the imaginary 
part of the response function S^^ojph, Qp/i, kj) which, 
in the case of one photon (u/p h, Qp/t), has been shown to 
read (see Eq. (16) of Ref. |38[ ) 



\n\ 2 L D 



E 



(44) 



|<r = 0,pi|77>| 5 



Jph 



2m e 



£(V,S=0) + 



(k,+Q ph ) 2 



2(2m e +m k ) 



where Q is the vacuum Rabi coupling, while is the 
relative motion momentum of the (e, X) pair. Since 
the momentum of the photocreated exciton is equal to 
the photon momentum Qp/j, the total momentum of 
the (e,A) pair is K; = k; + Q p h and its relative mo- 
tion momentum is such that k^ = + /3 e Ki and 
Qph = -Pi + /?xK 4 with /3 e = 1-/3 X = m e /(2m e + m h ); 
so pi = Px^-i — PeQph — /3xkj as the photon momen- 
tum is very small on the characteristic electron scale. By 
noting that 



2m e 2(2m e 



m h ) 



2(i 



rX 



(45) 



with fi e x being the effective mass of the electron-exciton 
pair defined above, we can rewrite Eq. (|4"4")l as 



SW(LJ ph ,ki)2!\Sl\ 2 L D 



|<r = O,0xki|»7}| s 



_ £(u,S=0) + Z 0+ 

(46) 



Before going further, we wish to make three comments: 
(i) In the absence of electrons in excess in the sample, 

the absorbed photon creates an exciton, the associated 

response function reading as 



sW(u ph ) = \n\ 2 L D J2' 



l(r = 0H| 5 



2(m c 



ih) 



i0+ 

(47) 



For bound exciton having a spatial extension ax , dimen- 
sional arguments lead to a normalized relative motion 
wave function (r|z/) such that 1 ~ o-x\( Y = ®\ v ) I 2 ! 80 



| (r = |z/>| 2 ~ (J- 
\ a x 



(48) 



It is worth noting that this quantity is of the order of the 
relative motion wave function for unbound trion states 
given in Eq. (|40"]) . This is quite reasonable because un- 
bound trions look very much like an exciton with a free 
electron moving possibly far away from it. 

We then note that Eq. gives an oscillator strength 
(clt I L) D smaller for bound trion than for unbound trion; 
for large sample size L, this a priori prevents drawing 
the full absorption spectrum for bound and unbound tri- 
ons using the same scale. It however is necessary to 
note that the L factor which here appears, comes from 
(u|p) = e ipu / 'L D < 2 ; so, this L factor physically is the 
length over which the particle keeps its momentum, this 
length being usually called "coherence length" . As a re- 
sult, when the sample volume increases, {ar/L) satu- 
rates to [clt I Lcoh) which can be not so small in poor 
samples. Nevertheless, we will here show the absorption 
spectra of bound and unbound trion states on different 
curves to avoid this scaling problem. 

(ii) The («0 + ) term in the response function of an exci- 
ton or a trion actually has to be broadened into ij with 7 
being associated with the inverse lifetime of the exciton 
or trion. Indeed, the response functions quoted in Eqs. 
(|44|47| follow from using the Fermi golden rule. This rule 
is known to be valid for transitions toward a continuum 
of states. When one photon Q p h is absorbed in an un- 
doped semiconductor, we form a well-defined exciton, its 
momentum being Q p h. So, the Fermi golden rule should 
not be used. We have shown (see Ref. [39.) that the 
Fermi golden rule can still be used for photoabsorption 
with exciton formation provided that the exciton level is 
broad enough to be seen as a continuum, the character- 
istic scale for this broadening 7 being the vacuum Rabi 
coupling Q. When the exciton level is very narrow, as 
in good quality samples, we say that the system suffers 
a "strong coupling" with the photon field: instead of ex- 
citons created by photon absorption through the Fermi 
golden rule, we have polaritons — the polaritons being the 
exact eigenstates of one photon coupled to one exciton. 
So, when using formulas like Eq. (|4"4"]) and Eq. (|4"T|) for 
photon absorption, we implicitly assume the trion or ex- 
citon to have a finite lifetime. As a result, i0 + should be 
replaced by ry. 

(iii) The linear response to a field having N p h photons 
is just Nph times the response function to one photon 
given in Eq. (|4"4")l . This makes the photon absorption lin- 
early increasing with photon number, i.e., laser intensity, 
as physically required for a linear response. 

Of course, in real experiments, we not only have more 
than one photon, but also more than one electron in the 
sample. It turns out that the effect of an electron number 
increase is far more subtle than a photon number increase 
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FIG. 13: (Color online) Absorption (oj p h) (defined in 
Eq. (|5ip ) in the presence of cold electrons in 2D quantum 

wells as a function of the photon energy uj p h in R^ D ^ unit. 
The trion ground state energy E^ ^ lies 0.24 R^ D ^ below a 



-4 R 



(3-D) 



(not shown). We 

2 , 



larger exciton peak which lies at 

have taken a mass ratio m/,/m e = 5 and set N p h\£l\* p = 1, 
the L D factor in Eq. (|51[) being canceled by the wave function 
part of this equation, due to Eq. (|39p . 



because photons are bosons so that they all have the same 
energy, while electrons are fcrmions; so, their energies 
spread out when their number increases. For sample with 
N e electrons, cold enough to possibly see trions, we must 
distinguish two regimes: 
(1) Cold carriers 

For temperature T small compared to the Fermi en- 
ergy ep for these N e electrons with ep — kp and N e ~ 
(kpL) D , the N e electrons are degenerate, i.e., they have a 
probability equal to 1 to have a momentum between 
and kp. So, the response function of these N e electrons to 
N p h photons should read, in the absence of interactions, 
as 



(49) 



with S^ T \ujph: ki) given in Eq. (|44j) . By replacing the k^ 
sum by an integral and by setting e = ((3 x ki) 2 /2p, eX = 
Pxkf /2m e , we get, for 2D systems having a constant 
density of states p, the photon absorption spectrum 



|(T) 



,UJph) 



-2 ImS^ T \co p h) associated with the trion 
ground state 770, through 

A%\u ph ) = 2nN ph \n\ 2 L D (50) 
pde (r = 0, y/2fx eX £\vo) ' 5(u> ph + e - £ {vo) ) 



fixe 1 



which readily gives, for a photon detuning 6 p h = £^ no ' 
uj p h positive, 



4oW) = 2nN ph \n\ 2 L D p (r = 0, ^2 t i eX S ph \ m ) 



for S^io) _ p X £ F < uj ph < £(»7o) an d zer0 otherwise. 

This shows that A^' {uj p h) essentially follows the de- 
crease of the bound trion relative motion wave function 
|(r = , p 1 770 ) 1 2 along p, the trion absorption spectrum 
spreading on the low energy side, since more energetic 
electrons can be used to form the ground state trion (see 
Fig. IS). 

It is worth noting that this absorption spectrum does 
not increase in amplitude under an electron number in- 
crease, but just spreads on the low energy side since kp 
increases when N e increases. Such a spectrum however 
is rather naive because it forgets many-body effects in 
the final state similar to the ones leading to Fermi edge 
singularities: Instead of having a smooth decreasing low- 
energy tail, these effects bring a peak at the absorption 
threshold, as discussed and observed in Ref. HfJ 

(2) Hot carriers 

In the case of hot carriers, the N e dependence of the 
absorption spectrum is totally different and strongly de- 
pends on where the carriers are injected. Let 



N(ki,e* e ) = N e 



E ki /(£k i)£ *)' 



(52) 



be their density probability. We will here use a Gaussian 
distribution 



/(£k i; e*) = e 



-C(kf/2m e 



(53) 



the constants £ and e* being possibly tuned through 
the bias voltage window of the injected electron current. 
Summing over the electron momentum kj then amounts 
to transforming Eq. (|4"9"]) into 

S^i^el) = Nph^NikusDS^iujp^ki). (54) 



1. Bound trion 

The absorption spectrum associated with the trion 
bound state 770, given by Eq. ([51]) in the case of cold 
carriers, transforms into 



Ag\uj ph ,e* e ) = 2nN ph \n\ 2 L D N e 



f(S ph /p x ,e* e ) 



(r = 0, \j2n eX 8ph\r)o) 



(55) 



(51) 



This shows that, for hot carriers, the absorption spec- 
trum is proportional to N e ; so, under an electron num- 
ber increase, the whole spectrum amplitude increases, 
instead of spreading out as for degenerate electrons. 

Figure [14] shows the 2D ground state trion absorption 
spectrum A^' {uj p h,£* e ) for various injected carrier aver- 
age energies e* (see Eq. ([53]) ). This spectrum can be 
used to probe the hot-carrier distribution. As the pho- 
ton energy decreases, the absorption intensity follows the 
decrease of the trion wave function |(r = , 1 77) | 2 with 
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FIG. 14: (Color online) Absorption Af^(u) p h,£t) defined in 
Eq. ()55|) associated with the trion ground state in 2D quan- 
tum wells, as a function of the photon energy Lj p h (in R^ D ^ 



unit) for different parameters el, and £ = 500 R^ D ^ in the 
electron distribution of Eq. (53) . We have set -AT pfl |S7| 2 iV e = 1 
and taken a mass ratio mk/m e = 5. 



relative motion momentum p;, as shown in Fig. 11(a) 



When the electron number increases, the amplitude of 
the whole spectrum increases, as a result of the N e fac- 
tor in Eq. (|55|) . but does not spread, in contrast to the 
low T limit shown in Fig. 1131 for which the absorption 
spreads without changing its amplitude. 



2. Unbound trion 

The situation for unbound trion is somewhat more 
complex because unbound trion states have energies close 
to the exciton energies; so, photoabsorption with forma- 
tion of a trion is mixed with photoabsorption with 
formation of an exciton A^ x > . The resulting absorption 
spectrum could be thought to read as 

fA^ + (1 - f)A™ = A^ + f [A^ - A<*>] , (56) 

where / is the capture rate of an electron by an exciton. 
For N e electrons in a sample volume L D , this capture 
rate should be of the order of the exciton volume divided 
by the average volume occupied by an electron, namely 



< =N ( a JL 
L D /N e e V L 



D 



(57) 



However, electrons with different kj's contribute differ- 
ently to the trion absorption as seen from Eq. (pT4"|) . So, 
the absorption spectrum should in fact read, instead of 
Eq. 421, 



A{u vh ,et) = + £iv(k ?; , £ 



k, 



x 'A^(cj ph ,ki)-A^\uj ph )\ (58) 



with N(ki,e*) given by Eq. (|52j) for hot carriers. We see 
that the absorption spectrum reduces to the exciton spec- 
trum in the absence of free carriers, N e = 0, as physically 
required. 

The exciton absorption spectrum A^ x '{uj p h) is made 
of delta peaks centered on the exciton energies e^">, and 
weighted by the value at r = of the exciton wave func- 
tion squared |(r = 0|z^)| 2 (see Eq. (g7J|). These delta 
peaks are broadened when taking into account the finite 
exciton lifetime, which amounts to replacing i0 + by 17 in 
Eq. (j47]) as previously explained; so, the delta functions 
are then replaced by Lorentzian function with a small 
half-width 7, namely 



8(u>) 



,2 + j2 



(59) 



Although, according to Eq. (|58|) . the photon absorp- 
tion contains exciton and trion parts, we have chosen 
here to only show the trion part in order to avoid am- 
biguity coming from their lifetime yet determined exper- 
imentally. The contribution to the photon absorption 
spectrum resulting from unbound trions is given by 

A^ {u ph , £ *) = iV p7l (^) ° ^ iV(k* , e* e )A^ (u ph , k, 



2TrN ph N e \fl\ 2 al 



x 



def(e//3 x ,e;) 
2 



Jo def{e/(3 x ,e* e ) 

x |(r = 0, y/2jl^\r)) ~ 6(u ph + e - £ M ) 



(60) 



the discrete sum over 77 being ultimately replaced by a 
continuous sum over k, ; . To compute it, we have intro- 
duced a finite lifetime and replaced the delta function by 
a Lorentzian function along Eq. (|5"9"|) . This essentially 
adds a small high-energy tail to the absorption spectrum 
which would have a sharp edge otherwise. 

Figures 15(a) and |15(bJ| show the lineshape of the three 
largest peaks located at the lowest s exciton levels for 

-4/9 R^°\ and 



(3D) 



2D quantum well, namely —4 R 

—4/25 Rx D \ these peaks being associated with unbound 
trions, i.e., the electron-exciton scattering states. We see 
that these peaks spread on both sides of each exciton 
level, due to energy conservation enforced by the delta 
function. Their peak height decreases with the increase 
of the exciton level due to the (r = 0\i/) factor in Eq. (|4"5|) 
which for 2D exciton is equal to 4/ ' ax\f2i:(2n — l) 3 with 
n = (1, 2, 3) for these three peaks. 

The peak lineshape essentially is a Lorentzian which 
moves slightly toward lower energy side when contribu- 
tions from electrons with large kinetic energy start to 
weight in. Moreover, the peak height slightly decreases; 
this can be attributed to broadening at large relative mo- 
tion momentum induced by electron-exciton scattering. 

Differences between calculated results and experimen- 
tal data are expected to primarily come from the fi- 
nite width of the quantum well. Indeed, all binding 
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twccn an electron and an exciton, we here construct a 
Schrodinger equation for trion using an electron-exciton 
basis, instead of the standard free carrier basis. One ob- 
vious difficulty with the electron-exciton formulation is 
to exactly handle electron exchange within the electron- 
exciton pair. This is done through exchange scatterings 
between electron and exciton similar to the ones which 
exist between two excitons, as appearing in the composite 
boson many-body theory. Restricting the exciton basis 
to the first three s-states, we have numerically solved the 
resulting Schrodinger equation for the trion ground state 
in 2D and 3D systems, as well as for excited bound and 
unbound states. Using just these three states, our re- 
sults for the ground state are in reasonable agreement 
with the ones obtained from the best variational meth- 
ods. The main advantage of the electron-exciton ap- 
proach presented here mostly is to reach excited states 
as easily as the ground state, these excited states being 
out of reach from usual variational methods. Since we do 
here solve a Schrodinger equation, we are also capable of 
reaching the trion wave functions of these bound and un- 
bound states. Through these trion wave functions, we 
can calculate the optical absorption spectrum associated 
with the trion bound states in the presence of injected 
hot carriers. The spectra have small peaks associated 
with the trion bound state and large peaks coming from 
unbound states, i.e., electron-exciton scattering states. 



VI. ACKNOWLEDGMENT 



FIG. 15: (Color online) Unbound trion absorption 
A^ T \ojph,e* e ) for 2D quantum well as defined in Eq. (|60[) 
where we have set N p h\Q\ 2 N e = 1. It shows unbound trion 
peaks resulting from electrons scattering with Is exciton (a) 
and with 2s or 3s exciton (b) as a function of the pho- 
ton energy uj p h in unit. We have taken a broadening 
7 = 0.01 R^ D ^ in Eq. (|59l) . Other parameters are the same 
as in Fig. [ft] 

energies are known to increase when the space dimen- 
sion D decreases — as evidenced from the fact that the 
exciton binding energy is equal to R^ D ^ in 3D but 

R-x^ = ^ m 2D. As a result, the trion binding 

energy in a finite-width quantum well is expected to be 
substantially smaller than its exact 2D value. Other siz- 
able effects are expected to come from the confinement 
of the trion wave function by potential barriers which 
are not infinite, and also by the potentials resulting from 
remote donor ions^ or quantum well interface thickness 
fluctuations^. 



V. CONCLUSION 

Motivated by the idea that Coulomb interaction be- 
tween an electron and a hole is far stronger than be- 
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Appendix A. SCHRODINGER EQUATION FOR 
TWO ELECTRONS AND ONE HOLE 

We here derive the Schrodinger equation © starting 
from Eqs. (|3l4p . The function 4>^'^) k appearing in this 

Schrodinger equation is linked to the V'kjf'k P re f ac t° rs 
appearing in Eqs. (|3I4I) through Eq. ((6j). The interest of 
this derivation is to evidence that the relevant function 
is not by no means V , k'if'k h but ^kl'k' kj,' 

We first consider the (S = 1,S Z = 2s) trion state given 
in Eq. ([3]). When used into the Schrodinger equation 

©, we get Eq. © with ^TkfL h replaced by ^kL' 
If we now project this equation onto (v\b Ph a p > : S >a Pe ;s, 
we find two terms that are no longer equal. However, 
since = U_ q , they can be combined to make appear 

(V'kTkfk^ " ^ZSj ' which is nothin g but 4Tkfk h for 
S = 1. The resulting equation then reduces to Eq. ©. 
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If we now consider the (S,S Z = 0) trion states given in 
Eq. we get the same equation ([5]) with a£ ^a^, , s 

replaced by [a] Ce . 1/2 a] ci ._ 1/2 - (-l) S aI c ;-i/2 a l<;i/2 
0ki5k replaced by . When projected 

onto (w|& Pfi a p ' ; _i/ 2 a Pc; i/ 2 , these two terms yield two 
terms that can again be combined to make appear 

(vi^k/kh + ( — ^) S ^kSkXh) ' 80 that wc a ^ so enc ^ U P w ith 
a Schrodinger equation for 0^ ^< ^ identical to the one 
given in Eq. ([9]). 



Appendix B. USEFUL RELATIONS ON 
ELECTRO N-EXCITON SCATTERINGS 



For completeness, this appendix contains the deriva- 
tion of a few useful relations involving excitons inter- 
acting with free electrons, some of which being possibly 
found in previous works on exciton^ and trion^ii. 

An exciton i = (V;, Qi) with center-of-mass momentum 
Qi, relative motion index Vi, electron spin s and hole spin 
m, expands on free electron-hole pairs as 



to an identity, as a result of the closure relation on \i) 
exciton states 



£a(*p-)A(p,T 



= Y, (^q^qMW'Ikeq.XqhPel*) 

i",Pe,qh,q' h 

= $i'iSk> e k c . (B.5) 

As seen from Fig. [2j for an exciton labeled by i = (t^, Q;) 
this Pauli scattering reduces to |p') (p|^i) provided 

that the (p,p') momenta are such that = p + a e Qi, 
k e = p' + a e Q;' and -p' + cthQi' = p + o^Q;- This 
gives 



A( .f e ) = 5k' < ,+Q i ,,k e +Q i (fi'|ke-aeQj')(k / e -a e Qj|fc'i). 

(B.6) 

These Pauli scatterings formally appear through two 
commutators 



B i';s',Bi 



D 



Si'iSs's - A's';is, (B.7) 

= ^E A ( 1? k ;>k;s- (B.8) 



t 



E «L A^keW 



k e k h 



E a p+"eQ,; S 6 -p+a h Q I ; m (Pl^>: (B.l) 



where a e = 1 — ah = m e j (m e + rrih). Conversely, a free 
electron-hole pair reads in terms of excitons as 



U k c ;s u k h ;m 



E^L^iw 



(B.2) 



E 



«,:.k. .k.:-.„, 



Vi\a h k e - a e k h ). 



From them, it is easy to show that the scalar product of 
two electron-exciton pairs is given by 



(vlap^Bj.^Bj ,al \v) 



(«|a P .; CT { 



B j;<?' > B i-s' 



B ts' B r^yi-M 



(B.9) 



Exciton and electron also interact by direct Coulomb 
scattering £ ( f f ) which also follows from two commu- 
tators. For an clcctron-holc Hamiltonian H , they appear 
through 



Since the hole spin m plays no role in problems dealing 
with one exciton interacting with an electron, we will, 
from now on, drop the index m from the notations. 

As easy to check from the above equations, electron 
exchange between one free electron k e with spin a and 
one exciton i with electron spin s leads to 



B l< 



= -EM^)4 ; .4 i;s > (b.3) 



k'i' 



where X( f kc ) is the Pauli scattering for electron ex- 
change shown in Fig. [TJ It readily follows from its di- 
agrammatic representation that this scattering is given 
by 



A( 



k' k e 



E)(i'|keP h )(p h k;|i) 



(B.4) 



This compact expression of the Pauli scattering between 
(e, X) pairs allows us to check that two exchanges reduce 



E^ X) B^ + V? 



E^K, 



(B.10) 
(B.H) 



Its diagrammatic representation, shown in Fig. [31 readily 
gives 



*(5 



k' k„ 



V K -k e E [^'l pe + ke - k^, P/; 



(i'\p e , Ph + k e - k' e ) ( PhPe |i).(B.12) 



For i = {vi,Qi), this scattering reduces to 
XX^i' Ip') (pWi) provided that the (p,p') momenta 
are such that p' + a e Qi> = p + a e Qi + k e — k' e and 
-p' + a h Qi< = p + Q'hQj. This gives 

€(^ k ;) = S Qil+KMi+K V K - K Y[(MP + *h(k e -K)) 



-(a h ->• -a e ) (p\ui) 



(B.13) 
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We can mix the direct Coulomb scattering £ with the 
electron exchange A to get the "in" and "out" exchange- 
Coulomb scatterings defined as 

= EWW?). (B.i4) 

k' e 'i" 

^° Ut (-^) = £f(? ^ )A( k :.' k I- (B.15) 

k'Ji" 

Using Fig. [51 we find that the "in" exchange scattering 
is given by 



k' k e 



[(i'|k e -q, P/l ) (B.16) 



q#0 PePh 

-(i'|ke q,Ph + q) 



If we now write the exciton i as (vi, Qj), Figure [6] gives 
this "in" exchange scattering as ^ Vq(i>i> |p') (p|i>i) pro- 
vided that the (p, p') momenta are such that k' e = 



p + a e Qi + q, p' + a e Qi> = k e - q, and -p' + a h Q ¥ = 
— p + ahQi for the electron-electron part, while we 
must have kg = p + a e Qi, p' + a e Qi' = k e — q, and 
— p' + a^Qi' = p + cn/iQ; + q for the electron-hole part. 
This ultimately leads to 



q 

(k' e - q- a e Qi|^) - (k' e - a e Q i |i/ i )l.(B.17) 



These two exchange-Coulomb scatterings are related to 
the dimcnsionless Pauli scattering for electron exchange 
through 



in / k' k e 



f out / k' k e 



(B.18) 

as easy to show by calculating (v\a-k' ;s Bi/. s HBl s al. . s \v) 
with H acting either on the left or on the right. 
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